Coherent Low-Energy Charge Transport in a Diffusive S-N-S Junction 
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We have studied the current voltage characteristics of diffusive mesoscopic Nb-Cu-Nb Joseph- 
son junctions with highly-transparent Nb-Cu interfaces. We consider the low-voltage and high- 
temperature regime eV < e c < ksT where e c is the Thouless energy. The observed excess current as 
well as the observed sub-harmonic Shapiro steps under microwave irradiation suggest the occurrence 
of low-energy coherent Multiple Andreev Reflection (MAR). 



PACS : 73.23.-b, 74.50.+r, 74.80.Fp 

A S-N-S junction made of a normal metal (N) embed- 
ded between two superconducting electrodes (S) exhibits 
a zero- voltage Josephson supercurrent [0. At the N-S 
interface, the microscopic mechanism is the Andreev re- 
flection of an electron into a hole which traces back al- 
most exactly the trajectory of the incident electron. This 
coherent process corresponds to the transfer of a Cooper 
pair in S and its inverse to the diffusion of an electron An- 
dreev pair in the normal metal. Let us consider a meso- 
scopic diffusive normal metal where the length L is larger 
than the elastic mean free path l e but smaller than the 
phase-coherence length L v . fn the long junction regime 
L ^> y^hD/A of interest here, the relevant energy scale 
is the Thouless energy e c = %/td = fiD/L 2 <C A. Here 
D = VFl e /3 is the diffusion coefficient, to = L 2 /D is the 
diffusion time, A is the superconducting gap. The essen- 
tial fact is that Andreev pairs with a low energy obey- 
ing e < e c remain coherent over the whole normal metal 
length L |2|]. This coherent window contributes signif- 
icantly to phase-coherent transport even if the thermal 
distribution width fc#T is much larger than the Thouless 
energy e c . This is exemplified in Andreev interferome- 
ters by the large magnetoresistance oscillations with an 
amplitude of about € c /ksT pj. Such a weak power-law 
temperature dependence is in striking contrast with the 
exponential decay of the Josephson coupling. 

At finite bias voltage (V), a dynamic regime takes 
place where the superconducting phase is time-dependent 
while successive Andreev-reflected electrons and holes 
gain the energy eV at each travel. This Multiple Andreev 
reflections (MAR) process leads to a sub-gap structure 
at the voltages 2A/pe (p integer) in the current-voltage 
characteristic |3|] and the energy distribution function [|4j . 
In quantum point contacts, the whole process of MAR 
is phase-coherent since the electron transit time is very 
short. This regime of coherent MAR was recently em- 
phasized in the description of both the current voltage 
characteristics, the shot noise and the d.c. supercurrent 
H of ballistic atomic contacts |i| |j. The related non- 
sinusoidal current-phase relation is indicative of coherent 
multiple charge transfer, but no direct observation has 
been reported so far. 



In a diffusive S-N-S junction, the diffusion of electrons 
through the normal metal involves the diffusion time td 
which must be compared to the period of the phase differ- 
ence x evolution Ty — h/2eV. The quasi-static condition 
for a small phase evolution during n successive Andreev 
reflections at each interface writes 2nrr) < Ty, which is 
equivalent to eV < Tre c /2n. In this regime, electrons and 
holes in the coherent window (e < e c ) can experience 
MAR while maintaining global coherence of the induced 
Andreev pairs. This process should imply the coherent 
transfer of multiple charges as in short ballistic point con- 
tacts, provided the interface transparency is close to 1. 
Recent experiments in diffusive S-N-S junctions focused 
on the opposite regime e c < eV ~ A pjl0||. Subgap 
structures were observed, but no signature of coherent 
MAR was found. 

The time evolution of the out-of-equilibrium energy 
distribution induced by MAR was also considered phe- 
nomenologically, taking into account the energy relax- 
ation. This lead to the prediction of an enhancement of 
the junction conductance at finite voltage |0,0|. The 
related a.c. Josephson coupling at twice the Josephson 
frequency was experimentally confirmed p3]. A close 
analogy can be drawn with the dissipative transport in 
a mesoscopic Aharonov-Bohm ring H] . Indeed both the 
relaxation current in such a ring and the MAR excess 
current in a S-N-S junction appear at bias voltages below 
h/r e 1 11 13J141 where r e is the energy relaxation time. 



In this Letter, we present an investigation of the co- 
herent dynamic regime in high-quality diffusive S-N-S 
junctions. We focus on the contribution of low-energy 
Andreev pairs to the current-voltage characteristic. We 
used as a probe an a.c. microwave field with a small pho- 
ton energy huj < e c so that it sits within the coherent 
energy window. Our main result is the observation of 
sub-harmonic Shapiro steps that we discuss in terms of 
coherent multiple pair transfer at low energy. Here, the 
relevant number of Andreev reflections Hto/2eV is im- 
posed by the photon microwave energy, not by any energy 
gap. 

The samples consist of a small Cu conductor attached 
to two Nb electrodes which are superconducting below 



T c ~ 8 K (Fig. |l|a inset). The electrical contacts in- 
clude an on-chip capacitance (0.2 pF) which connects 
the microwave circuit made of a cryogenic 50 51 coaxial 
cable to the sample. The fabrication process makes use 
of e-beam lithography and a shadow evaporation tech- 
nique based upon a thermostable resist compatible with 
UHV electron-beam Nb evaporation |15| ], The two sam- 
ples described here (labelled A and B) belong to the se- 
ries of S-N-S junctions investigated in Ref. jlq ]. They 
show a similar behavior. We intentionally kept the ther- 
mal fluctuations negligible at every temperature by de- 
signing wide junctions so that their normal-state conduc- 
tance Gn and hence the Josephson energy Ej = hl c /2e 
is large : Ej > k,BT c . The sample and measure- 
ment parameters were also chosen to fulfill the condition 
eV < e c < ksT < A. The parameters for sample A 
(respectively B) are the following : length of the normal 
metal L = 710 (800) nm, width 580 nm and thickness 100 
nm. The diffusion coefficient is D = 250 (230) cm 2 /s and 
the normal state resistance is 0.152 (0.183) fl. The cal- 
culated Thouless energy e c = 33 (24) /j,eV coincides with 
the values obtained from the magnitude and tempera- 
ture dependence of the d.c. critical current I c assuming 
perfect transmission at the S-N interfaces Jig] . 

Fig. pi shows the differential resistance of sample A 
as a function of the d.c. current bias. The critical cur- 
rent I c is easily identified by the sharp jump in differ- 
ential resistance, even at high temperatures when it is 
strongly reduced. The differential resistance character- 
istics shows a striking behavior above 4 K : the ampli- 
tude of the main peak at the critical current decreases 
while a broader bump develops at higher current. In 
this temperature range the current-voltage is fully re- 
versible, which discards heating effects. The comparison 
in Fig. |ljb with the ohmic behaviour and the resistively 
shunted junction (RSJ) model V — Gjj ^fP — T* (at 
5.5 K) demonstrates that this behavior corresponds to 
a low-energy excess current. The observed features are 
reminiscent of the foot-like structure previously observed 
in superconducting micro-bridges [|i7| , |l8[ which was as- 
signed to non-equilibrium processes in presence of a gap 
oscillation. In our long S-N-S junctions, the role of the 
gap is played by the mini-gap which sets the energy scale 
at e c . The microscopic origin of the excess current is the 
coherent MAR : each Andreev reflection at the N-S in- 
terfaces transfers a Cooper pair in and out, resulting in a 
current increase until an inelastic process takes place [|19[ . 
In average, each coherent electron transfers T e /To times 
the elementary charge. From the non-equilibrium model 
of Ref. 113] , a differential resistance bump is predicted at 
eV = h/r e . In the experiment, the position of the dif- 
ferential resistance maximum corresponds to r e — 90 ps 
at T = 5.5 K. This value is consistent with the known 
electron-phonon relaxation time in Cu pO| and exceeds 
the diffusion time tjj — 20 ps at every temperature. It 
gives a phase-coherence length L v « 1.5 /im ]2l|] which 



is about twice the normal metal length. The amplitude 
of the conductance enhancement dG was estimated by 
Zhou and Spivak Q as : 3G/Gn — ec/ksT . t e /t£i. 
This gives here 8G/Gm — 0.3 at 5.5 K, in qualitative 
agreement with the experiment. 
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FIG. 1. a) Measured differential resistance of sample A as 
a function of d.c. bias current at different temperatures. Inset 
: Electron micrograph of a typical sample made by shadow 
evaporation, b) Current-voltage characteristics as obtained 
by numerical integration of differential resistance curves. Also 
shown is the theoretical curve calculated from the RSJ model 
at 5.5 K (dotted line) and the ohmic behaviour (dashed line). 
The data all sit in the coherent window eV < e c = 33 /xeV. 

We investigated the coherent nature of the low voltage 
electrical transport by applying a low power (-20 dBm) 
microwave field. This low power ensures us that the mi- 
crowave acts only as a probe and does not drive the elec- 
tron energy distribution. Fig. shows the differential 
resistance and the current as a function of voltage for 
sample B at 4 K in the presence of an a.c. current I u 
of frequency 6 GHz. The feature corresponding to the 
Shapiro step is observed at voltage V\ — TiLu/2e = 12.4 
//V (index 1 on the figure), as expected. In addition, we 
observe two structures at one half : Vl/2= 6.2 //V and 
one third : V1/3 = 4.1 /iV of this voltage (index i and 
■|). The latter is clearly visible only on the differential re- 
sistance curve. Because of the high temperature and the 
small a.c. excitation, the plateau amplitudes are small 
and rounded by thermal fluctuations. The magnitude of 
these structures was observed to increase linearly with 
the microwave current I u , leading to true plateaux with 



zero differential resistance at high a.c. power (-10 dBm). 
Then, a large number of peak structures, both harmonic 
and sub-harmonic, emerge from the noise level. We found 
that their location is independent on both the tempera- 
ture and the microwave power. The position of the ob- 
served peak structures as a function of the microwave 
frequency between 4 and 18 GHz is plotted in Fig. ||. 
The harmonic 1 represents the fundamental Shapiro step. 
The ensemble of the peaks obeys the simple law given by 
V = m/n . fru>/2e. While the integer steps (n = 1) are 
simple results of the ordinary Josephson coupling, the ob- 
servation of fractional steps (n > 1) deserves a detailed 
discussion. 
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FIG. 2. Lower curve (right scale) : Differential resistance 
of sample B at 4 K as a function of voltage in presence of a 
microwave current of frequency 6 GHz and power -20 dBm. 
Structures induced by the microwave current appear at the 
usual Shapiro step position Vi = hui/2e = 12.4 juV and also 
at one half and one third of this value. Upper curve : the 
corresponding voltage steps on the current-voltage character- 
istic. 

In order to extract the step amplitudes, we assume that 
the steps can be described by the RSJ model : the step 



width is twice an effective critical current 



i/i 



smoothed 



by a thermal broadening. We used either an analytical 
expression or the high temperature approximation 0] to 
determine for each temperature and each step the actual 
half width i m / n . Fig. |J summarizes its temperature de- 
pendence for the three main Shapiro steps (1, i and ^) 
at 6 GHz microwave frequency. At the center of the step 
1 plateau, the current bias / ~ Gn V\ is much larger than 
the critical current I c , so that the step amplitude should 
follow the usual voltage-bias law i\ = I c Ji{Iu/I) [p2| - 
Here J\ is the Bessel function of the first kind. As can 
be seen in Fig. ^ the temperature dependence of the am- 
plitude i\ strictly follows that of I c . The small observed 
ratio i\/ I c ~ 0.06 allows an estimate of the a.c. excitation 
current I u w 10 /iA at -20 dBm. The treatment of the 
microwave current as a non-perturbative probe is hence 
fully justified. The most striking observation is the non- 
monotonic temperature dependence of both ii/ 2 and ii/3, 
in strong contrast with the monotonous behavior of i\. 



As seen in Fig. [|, i 1 / 2 starts increasing at temperature 
of about 3 K, then reaches the same amplitude as i\ near 
4 K, and finally slowly decreases at higher temperatures. 
Let us point out that the temperature at the maximum 
is much larger than the Thouless temperature e c /fcs = 
0.28 K. A similar behavior is found on the smaller ^ step. 
The persistence of the ^ and 3 at high temperatures, i.e. 
when the d.c. Josephson coupling is suppressed, is the 
most important result of this work. 
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FIG. 3. Peak position as a function of frequency for integer 
(black symbols) and non-integer steps (open symbols) in sam- 
ple B. These steps were observed at a 4K temperature and a 
-10 dBm microwave power. 

The persistence of (integer) Shapiro-like steps at high 
temperatures |23J is expected to occur in a four-terminal 
configuration under current injection from a normal 
reservoir, but not in our two-terminal geometry. One 
can also discard the inductance or the capacitance of the 
junction itself or the onset of chaos as a possible ori- 
gin of sub- harmonic steps pj, as our high- conductance 
junctions are strongly overdamped. The sub-harmonic 
steps must indeed be viewed as a manifestation of a non- 
sinusoidal Josephson current-phase I(x) relation. An in- 
trinsic non-sinusoidal contribution to the d.c. supercur- 
rent is expected at zero temperature but strongly sup- 
pressed as the temperature increases J24J . In a long junc- 
tion at ksT > 5 e c (k2K for our samples) the second 
harmonic is expected to contribute to less than 2 % of 
the total current. In contrast, the sub- harmonic steps ob- 
served here not only appear as the temperature increases, 
T > 3 K as seen on Fig. [|, but they eventually exceed 
the integer steps. The intrinsic non-sinusoidal phase de- 
pendence of the d.c. Josephson current therefore cannot 
explain the sub-harmonic steps. The non-equilibrium a.c. 
Josephson coupling model ||12|| introduced in our discus- 
sion of the excess current was shown to predict a finite 
sin 2x term even in absence of d.c. Josephson critical cur- 
rent. Lehnert et al. p3[ reported such observation in InAs 
based junctions in the clean limit. Our results in the dif- 
fusive regime show strong similarities with fll3[ ] although 
the latter ones are concerned with higher electron energy 



scale because of the smaller electron transit time. 
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namics of superconducting nanocircuits" . We thank F. 
Hekking for pointing out the relevance of coherent MAR. 
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FIG. 4. Temperature dependence of the half width i m / n 
of integer and fractional steps in sample B under -20 dBm 
microwave irradiation at 6 GHz. The detection threshold was 
about 0.1 fiA. On the same graph is plotted the measured d.c. 
critical current 7 C . 

Let us now discuss the implications of sub-harmonic 
steps in terms of coherent MAR processes. The bias volt- 
age involved here is so low that the phase difference is 
quasi-static and low-energy Andreev electron pairs are 
phase-coherent over the junction length. This is the 
regime where electrons close to the Fermi level experience 
coherent MAR. The | step is therefore associated with 
the coherent transfer across the diffusive metal of double 
pairs. Along the same way, the i step corresponds to the 
transfer of triple pairs. In the case of an interface trans- 
parency close to 1, this interpretation holds provided the 
phase-breaking events can be neglected. In this respect, 
the characteristic time scale 2nTjj = 2nx 28 ps has to 
be smaller or at least comparable to the phase coherence 



time t v ~ 160 ps at 4 K |20|. At this temperature, the 



condition is fulfilled for n = 2 and 3, which is consistent 
with our observations of ^ an d 3 steps at low power. 
Higher-order steps up to n = 5 arise only at high power 
when the inelastic damping is compensated by a larger 
a.c. current. 

In summary we have investigated the low-voltage co- 
herent transport in diffusive S-N-S junctions with high 
barrier transparency. We observed a series of sub- 
harmonic Shapiro steps at high temperature and low volt- 
age eV < e c < ksT where ordinary Josephson coupling 
is known to vanish. Presently there is no satisfactory 
description of the microscopic mechanism which is very 
likely to involve coherent transfer of multiple charges. 
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